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ABSTRACT 



Two experimental techniques to determine the frequency- 
dependent complex mass of a viscous fluid contained in a rigid porous 
solid are investigated. In one technique the moving mass of the fluid 
is sensed by its effect on the measured input electrical impedance of a 
moving coil transducer. In the second technique the moving mass is 
extracted from the measured acoustic pressure required for the fluid 
to oscillate with a known amplitude through the solid frame. A 
description of the apparatus and preliminary results using the 
pressure method as well as a comparison of the two techniques are 
presented. 
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I. INTRODUCTION 



A. BACKGROUND 

Acoustic propagation in the ocean is frequently controlled by the 
sea bottom (Urick, 1982, p. 11-1). Reflective bottom loss due to 
impedance mismatch, scattering due to bottom roughness, and 
attenuation due to dilitational and shear wave excitation in bottom 
sediment all dominate shallow water transmission at low frequencies, 
and influence deep water propagation to short and intermediate 
ranges at mid-frequencies (Tolstoy and Clay, 1972, p. 8). The 
importance of these physical phenomena to underwater acoustics is 
well recognized. Extensive literature is available on the broad subject 
of marine sediment acoustics because of the vast contributions by 
many renowned individuals over the past several decades, e.g., 
Hampton (1974), Hamilton (1974), Kuperman and Jensen (1980), 
Akal and Berkson (1986), and Stoll (1986). 

Even so, one example of bottom interaction physics that remains 
difficult to specify theoretically and validate experimentally is the 
frequency dependence of the attenuation. Because of the difficulty of 
measuring attenuation at low frequencies, values are obtained by 
extrapolation from accurate high frequency measurements; too often 
decades removed from the frequency span of study (Hovem. 1980, p. 
2). This problem, coupled with the inherent experimental 
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uncertainties and the fact that no single data set spans the entire 
frequency range of interest in underwater acoustics, i.e., Hz to MHz 
(Kibblewhite, 1989, p. 718), underscores the requirement for an 
accurate theoretical model for the attenuation of sound in sediments. 
Such a theoretical model can be prescribed by applying the Biot theory 
of acoustic wave propagation. Biot theory accounts for fluid-solid 
viscous interaction, which is an important ingredient in the 
phenomena of bottom attenuation. However, the details of this 
interaction are not well characterized for porous solid geometries. 

The Biot model of fluid-saturated sediment contains thirteen 
parameters, defined in Biot (1956a, 1956b), Stoll (1974,1985,1986), 
and Plona and Johnson (1984), which depend on the microscopic 
pore geometry of the porous solid, the bulk properties of the solid, 
and the bulk properties of the fluid. These parameters are very 
difficult to extract from conventional measurements of compressional 
sound speed and attenuation. This fact, along with the paucity of 
independent observations of the effects of fluid-solid interaction, has 
been one of the major obstacles to general acceptance of the Biot-Stoll 
view by the underwater acoustics community (Kibblewhite, 1989, p. 
718). However, recent developments include success in defining 
some of the key Biot parameters and coefficients. As reported by 
Kibblewhite (1989); Bedford et al (1984), Johnson et al (1987), and 
Baker (1986) have contributed significantly to specifying parameters 
required as input to the Biot formulation of elastic wave propagation in 
a fluid saturated porous solid. 
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Least well known are the three parameters which depend on the 
microscopic pore geometry and which, along with the bulk fluid 
density, completely determine the so-called complex mass density, P , 
of the fluid in the pores (Ogushwitz, 1985, p.430, and Stoll, 1985, 
p. 1798)(as revealed in Chapter II, the complex mass density for a 
rigid and stationary porous solid is the negative of the ratio of the 
gradient in fluid pressure to its acceleration). For this reason, a 
measurement of P is a superior probe of these parameters (Baker, 
1986, p. 2). 



B. OBJECTIVES 

The goal of this research program is to determine the frequency- 
dependent complex mass of a viscous fluid in a rigid and stationary 
porous solid for various typical pore geometries. Two experimental 
techniques are being investigated, with the flu id -saturated solid 
sample contained in a small cylindrical cavity, depicted in Fig. 1 
(Grant, 1988, p. 4). The first technique, the so-called impedance 
method, is based on electroacoustic network theory. The moving 
mass of the fluid is sensed by the measured input electrical impedance 
of moving coil drivers. Preliminary results for the impedance method 
have been reported by LCDR Grant (Grant, 1988, p. 57). The second 
technique, the so-called direct method, is based on Newton’s Second 
Law for a viscous fluid contained in a rigid porous solid. Stiff 
piezoceramic drivers are used to apply a known displacement for a 
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given voltage input. The moving mass is extracted from the observed 
pressure difference across the sample required for the pore fluid to 
oscillate with a known amplitude. 
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Figure 1. Simple Model. 

The specific objectives of the research described in this thesis are: 
(1) to develop and calibrate an apparatus to extract the complex fluid 
mass contained in a rigid porous solid from acoustic pressure and 
particle velocity measurements (direct method), (2) to compare these 
calibration results with those of the impedance method (Grant, 1988, 
pp. 44-56), and (3) to identify the better method of the two for 
further investigation with various porous solid sample geometries. 
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II. THEORY 



A. BIOT-STOLL THEORY OF SOUND PROPAGATION IN A FLUID- 

FILLED POROUS SOLID 

As summarized by Stoll (1974), propagation of sound in the ocean 
bottom is mainly controlled by two fundamentally different 
mechanisms: the first is the friction at grain-grain contacts due to the 
inelasticity of the solid frame, and the second is the viscous drag 
between the pore fluid and solid frame. Sediments of low 
permeability also exhibit another kind of viscous damping, attributed 
to local fluid motion around intergranular contacts (Stoll, 1985, p. 
1790). 

In principle, the motion of solid and fluid components can be 
described by differential equations for a microscopic volume element 
of porous material, subject to the appropriate boundary conditions at 
the solid frame-pore fluid interface. Solutions to these equations can 
be found only if the microscopic pore geometry is completely 
specified, a practical impossibility. However, a phenomenological, 
macroscopic (averaged over the microscopic structure) description of 
motion is possible, as long as the acoustic wavelength is much greater 
than a characteristic pore radius. Biot derived just such a description 
of motion in his series of classic papers (Biot, 1956a, 1956b, 1962). 
Many equivalent formulations of Biot equations and parameters exist to 
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suit particular applications. My choice is the set of equations of 
motion presented by Baker (1986, p. 29): 



Pi^-u = -p a ^(u -U) + B V(V u) + c V(V U) 
dt 2 dt 2 

-NVxVxu - bF^lu - U) 

^ ( 2 . 1 ) 



p 2 |i U = -p Au - u) + cv(v- u ) + DV(V- U ) - bF|(U - u) 
dt 2 dt 2 * (2.2) 



where, 



U = average fluid displacement, 
u = average frame displacement, 



P l = mass density of bulk solid multiplied by the volume 
fraction of solid, i.e.. Pi =(1 - P)p s , 

p 2 = mass density of bulk fluid multiplied by the volume 
fraction of fluid, i.e., P 2 = Ppf , 

p a = effective mass density normalized to the apparent 
fluid density, i.e., Pa = (a - l)Ppf , 

P = porosity, defined as the volume fraction of fluid, 

a = tortuosity parameter, a function of the microscopic 
pore geometry, 

b = friction coefficient = PRdc '» Rdc Is the DC flow 

resistance, defined as the average gradient in fluid 
pressure divided by the average fluid velocity in the 

direction of -^P, 

N = shear modulus of the frame, 
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F = complex, frequency-dependent, universal correction 
function. 



The coefficients B, C, and D are generalized elastic constants, and are 
functions of the bulk fluid and solid elastic contants, and the porosity, 
as follows: 



B = 



l-P-^i-PJK. + P^Kb 



lp i7. 



+ i k7 



+ t n 
3 



(2.3) 



where, 





PK S 


l -p 


+ P^ s - 


K S J 


K f 


p 2 k s 


1 -P-^. 


+P *k 


K S J 


Kf 



(2.4) 



(2.5) 



K s = solid grain bulk modulus, 

K b = solid frame (absent the pore fluid) bulk modulus, 



Kf = pore fluid bulk modulus. 

For the case of a stationary frame, to which this study is restricted, 
u = 0. Noting that pa + p 2 = «Ppf , Eq. (2.2) yields, for u = 0: 



ocPpr^-U = DV(V • U ) - bF^ 

at 2 at . (2.6) 

For a perfectly rigid frame Kf « K s ,Kb t and Eq. (2.5) reduces to 
D = PK f . Recalling also that b = PRdc . Eq. (2.6) becomes: 
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(2.7) 



aPpt^-U = PK f V( V • U ) - PRdcF- ^ 

Assuming harmonic time dependence (ei® 01 ) , Eq. (2.7) becomes: 

-to^apf + ^f|U = K f V(v • U ) 

From the linearized continuity equation, 



,dU 



( 2 . 8 ) 



VU =-?£l 
Pf 



(2.9) 



Applying this result and the definition 



K ( = Pf SL 

5p f 



for the bulk modulus to Eq. (2.8) gives: 



. w 2| a p f + ^ F j U = _Vp 



Defining the complex mass density by 



p = apf + 

J« 



( 2 . 10 ) 



( 2 . 11 ) 



( 2 . 12 ) 



Eq. (2.9) can be cast in terms of fluid particle velocity u=j<oU as 

jcopu = -Vp ( (2.13) 

which is in the familiar form of Euler's equation (KFCS, 1982, p. 104), 
(hereafter the subscript f on P will be omitted; it is to be understood). 
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The real part of Eq. (2.12) is the apparent fluid mass density and 
the imaginary part is the negative of the apparent flow resistance 
divided by the angular frequency: 



Both Eq. (2.14) and Eq. (2.15) depend on frequency. The normalized 
complex mass density, P / Pf , is called the complex, or dynamic, 
tortuosity, a . and depends only on the microscopic pore geometry 
and frequency. In the context of the Biot theory, the microscopic 
pore geometry is characterized by the three parameters: (intrinsic) 
tortuosity a, DC flow resistance Rpc. and the structure factor 5, which 
is a scale factor for the argument of the complex universal correction 
function F. In Biot's original notation, these correspond to the pore 
size parameter, permeability, and structure factor, respectively 
(Ogushwitz, 1985, p. 430). 

Viscosity of the pore fluid is very important to solid/fluid coupling. 
A pertinent quantity, therefore, is the viscous penetration depth (5 V ), 
or skin depth, defined as 



Re(p) = p ef f{to) = Re| 




(2.14) 




(2.15) 




(2.16) 
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o 

where rj is fluid (shear) viscosity, (e.g., 5 V = 56p, for water at 20 C and 
to = 100 Hz). The relative magnitude of the imaginary and real parts of 
the complex mass density, and therefore the relative importance of 
fluid viscosity versus inertia in determining the fluid motion, depends 
upon the ratio of the viscous penetration depth to a characteristic 
pore size. If 8 V is small compared to a characteristic pore size, then 
the motion of the fluid is very similar to that of an inviscid fluid, and 
independent (of the solid) sound propagation within the fluid is 
possible. This case is termed the high frequency limit. By contrast, if 
the viscous penetration depth is large compared to the characteristic 
pore size, the fluid is essentially locked to the rigid solid frame. In 
this case, sound cannot propagate independently within the fluid, and 
it is referred to as the low frequency limit. (Bourbie et al, 1987, p. 
86 ) 

Two methods for experimentally determining the complex mass 
density of a flu id- saturated porous solid are being investigated, termed 
the direct method and the impedance method. A theoretical 
description of each follows. 



B. DIRECT METHOD 

In the so-called direct, or Newton's Second Law, method, the 
complex mass density p is directly obtained from a measurement of 
the acoustic pressure gradient required to produce a certain (known) 
acceleration, using Eq. (2.13): 
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p = 



:Ye 

jwu . 



(2.17) 



The description of how pressure and acceleration are measured is 
given in the next chapter. 



C. IMPEDANCE METHOD 

In the impedance method, the complex mass density is indirectly 
obtained from a measurement of the input electrical impedance (Zm) of 
a moving coil transducer as follows, from Wilson (1985, p. 28): 

zin=ze+ 2i£ _ 

R m+ R t + j<o(M + H) + j t (218) 

where; 

B1 = transduction coefficient of the driver, 

R m = mechanical resistance of the driver, 

R r = effective mechanical resistance of the fluid, 

M = moving mass of the driver, 

Mr = effective mass of the fluid, 
s = stiffness of the driver, 

Ze = blocked electrical impedance of the driver, 

The effective mass, M r , and the effective mechanical resistance, R r ,of 
the fluid, obtained from the imaginary and real parts of Eq. (2.18), 
respectively, are given by: 



1 1 



(2.19) 



H = £lm 



(Bl) 2 
Zin - Ze 



+ - a - - M 

CO 2 



and 



R r = Re 



(Bl l ' 
Zin - Ze. 



( 2 . 20 ) 



Finally, Mr and R r are related to P eff and R fiow , respectively, by 



_ My p _ R r 

Peff PSL and flow PSL , (2.21) 

where; 

P = the porosity, 

S = the cross -sectioned area of the chamber, 

L = the length of the fluid plug, 

PSL = the volume occupied by the fluid. 



A description of the apparatus developed to extract the complex mass 
from impedance measurements using a moving coil loudspeaker 
follows in the next chapter. 
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HI. EXPERIMENTAL DESIGN. PROCEDURE. AND PRELIMINARY 

RESULTS 



A. BACKGROUND AND PREVIOUS WORK 

In Chapter II, theoretical relationships were derived between 
physically measurable quantities by which the complex moving mass of 
a fluid contained in a rigid porous solid can be extracted. This 
chapter describes the design and preliminary results of the two 
experimental techniques being investigated for this purpose. 

For his Master’s research, LCDR Grant developed the apparatus for 
the impedance method, (Grant, 1988), and his work is summarized 
here. The development of the impedance method was accomplished 
in two phases: (1) selection and characterization of two matched 
moving-coil transducers, and (2) design and construction of a test 
chamber. A pair of Philips model AD062T8 one-inch diameter 
tweeters, depicted in cross-section in Fig. 2 (Grant, 1988, p. 18), 
were used as drivers. To ensure oscillatory plug flow of the fluid 
through the sample of porous material, the pair of transducers were 
driven in parallel and out-of-phase with respect to each other in a 
system closed to ambient pressure. To satisfy the rigid frame 
condition, a massive six inch outer diameter brass test chamber, with 
a centered 1.5 inch diameter cylindrical cavity, was designed and built 
as shown in Fig. 3 (Grant, 1988, p. 40). Porous samples are to be 
housed within a cylindrical metal sleeve, and the sleeve inserted into 
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Voice Coil 

Figure 2. Philips Driver. 



the cavity. Overall test chamber length was limited to three inches to 
avoid acoustic resonances in the experimental frequency range of 
interest, i.e., below one kilohertz. 



Vent 

Connection 



BMC 




Gas 

Connection 



Figure 3. Cutaway View of Impedance Method Apparatus. 
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In a series of experiments, with the objective to calibrate the 
apparatus with the sample cavity filled soley with water, LCDR Grant 
successfully demonstrated that the moving mass load presented to the 
face of one moving coil transducer could be extracted from its 
measured input electrical impedance. This required a meticulous and 
thorough investigation of the five properties of the moving coil driver 
which appear in Eq. (2.18) (the details of which are presented in 
Grant, 1988). Reported experimental uncertainty in all five quantities 
was consistently less than 4.5 per cent (Grant, 1988, pp. 21-33). 
High-precision input electrical impedance measurements were made 
using a Hewlett-Packard Low Frequency Impedance Analyzer Model 
4192 (Hewlett-Packard Corporation, San Jose, Ca.). LCDR Grant 
concluded that it appeared possible to extract the effective mass with 
an accuracy of about five percent over a significant portion of the 
frequency band 100 to 1000 Hz; e.g., Fig. 4 (Grant, 1988, unpublished 
results). 
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Radiation Mass for 20 ml of Added Water 




Figure 4. Impedance Method Mass Extraction for 20 ml of Added Water. 



B. FURTHER IMPEDANCE METHOD EXPERIMENTS 

LCDR Grant's results using the impedance method revealed a 
systematic error with increasing frequency in the band of interest, as 
Fig. 4 indicates. The magnitude of this error was found to be greater 
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for larger added mass. A considerable effort was therefore made in 
the present investigation to improve the impedance technique. A 
discussion of the source of the error and a description of the efforts 
made to correct it follow. 

Fig. 5 (Grant, 1988, p. 15) and Fig. 6 (corrected figure from Grant, 
1988, p. 14) show the block diagram and lumped element 
representation of the equivalent electrical circuit for the impedance 
method. The subscripts M and R indicate the electrical equivalents of 
the transducer open circuit mechanical impedance and the radiation 
impedance, respectively, in the mobility analog (KFCS, 1982, p. 360). 
In this analogy, a greater mechanical or radiation impedance results in 
a smellier equivalent electrical impedance. For frequencies below one 
kilohertz, the impedance of the inductance (representing mechanical 
compliance) is greater than any other element, and this component is 
largely ignorable. Hence, for increasing frequency, the radiation and 
mechanical impedances appear more and more as short circuits. 
Under these conditions, the degree to which the blocked electrical 
impedance is known determines the precision with which the 
radiation impedance can be extracted from the measured input 
electrical impedance. 

LCDR Grant realized the blocked electrical impedance cannot be 
represented as a simple resistor and inductor in series (as in Fig. 6). 
In fact, a frequency dependent blocked electrical impedance was used 
to obtain the results shown in Fig. 4. Because it is virtually impossible 
to completely block a transducer without damaging it, however, the 
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value of the blocked electrical impedance which Grant used in his 
analysis was obtained from electrical impedance measurements on a 




Figure 5. Equivalent Circuit. 
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Figure 6. Lumped Element Equivalent Circuit - Mobility Analogy. 



similar moving coil transducer which had been completely potted with 
epoxy (Grant, 1988, p. 22). For the extraction of moving masses of 50 
g or so, the blocked electrical impedance must be known with much 
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more certainty, i.e., for the particular driver being used. 
Furthermore, because the mechanical resonance of the driver with no 
added mass occurs at about 1000 Hz, some of the error in the 
extracted value for the larger masses could be due to insufficient 
knowledge of the mechanical impedance. In summary, it is believed 
that extension of the impedance method to successfully extract 
masses of 50 g or more hinges on knowing the transduction 
coefficient, mechanical resistance, moving mass of the driver, stiffness 
of the driver, and the blocked electrical impedance, to a higher 
degree of accuracy. 

The first effort to improve the impedance method involved finding 
the best values of these parameters, that is, the values which gave the 
best results for the extracted mass for a collection of known added 
masses. A computer program was written to control an HP 4192 to 
measure the input electrical impedance. Impedance measurements 
were made for added masses of 0 to 100 g in 10 g increments from 10 
to 2000 Hz in 10 Hz steps. A grid-search least squares program 
adapted from Bevington (1969, pp. 208-213) was used to find values 
of the "known" transducer parameters (Bl, Zm. Ze) which minimized 
the magnitude of the fractional error in the computed and actual 
radiation impedance (jGJMd, 



(Blf 

(Zin - Ze) 






jcoM r 



2 



1 



(3.1) 
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and to store the results In a file for use by graphics routines. Seed 
values for the parameter values were taken from Grant (1988, pp. 
13,39). The results showed no significant improvement over Grant's 
results. 

As an alternative attempt to improve the impedance method, the 
problem of finding the optimum values of the transducer parameters 
in Eq. (2.18) was re-cast so as to find the complex, frequency 
dependent transduction coefficient, Bl, that gives the least variance of 
the mechanical impedance of the transducer extracted from 
impedance measurements. From Eq. (2.18), the mechanical 
impedance may be obtained from measurements of input electrical 
impedance with a mass added, Zin , and with zero mass added, 
^in, no mass , according to: 



Z™ = 



• IN, NO MASS 



- Zin 



(Bl) ' 



(Zm + jcoM r ) 



-1 



which has the root 



(3.2) 



T 7 _ jwM r 
— ~ t 

2 



(jcoM r ) 






4 (Z IN, NO MASS - Zin)J . (3.3) 

This second attempt to obtain better values for the transducer 
parameters produced results within the experimental error reported 
by Grant (1988, p. 30) in the band 500 to 1500 Hz. For frequencies 
greater than 1500 Hz, however, the optimum value of Bl was 
unreasonably low, for reasons unknown. Following the disappointing 
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results of attempts to improve the impedance method using the 
available apparatus, it was concluded that better results were not likely 
to be obtained without changing to another driver. Further 
investigation focused on the development of the direct method, 
described next. 



C. DIRECT METHOD DESIGN 
1. Closed System Design 

Figure 7 shows a cutaway view of the new apparatus designed 
and built to extract the complex fluid mass from acoustic pressure and 
velocity measurements (direct method). The sample chamber is 
unchanged from the impedance method apparatus. The direct 
method requires stiff transducers, driven in such a way as to cause the 
fluid in the sample to oscillate as a plug. For this reason, a "stack’’ of 
piezoceramic rings was chosen for the motor element of each driver. 
Axially polarized Navy Type I lead zirconate-titanate rings, .75 inches 
in diameter and .25 inches thick, with a .25 inch diameter center 
hole, were procured (Channel Industries, Inc., Santa Barbara, Ca.). 
The piezoelectric strain constant in the direction of polarization, CI 33 , 
was measured for each element using a Berlincourt CI 33 meter (Channel 
Products, Inc., Chesterland, Oh.). From four selected ^-matched 
pairs, two nearly identical piezoceramic transducer stacks were 
assembled by cementing together four rings for each driver. 
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A piston and surrounding diaphragm assembly seals the cavity at 
either end of the sample. An impedance head, B&K Model 8001 
(Bruel & Kjaer, Foster City, Ca.), was positioned between the 
piezoceramic stack and this assembly. An impedance head is a 

Vent Piezoceramic "stack" 




Figure 7. Cutaway View of Direct Method Apparatus. 

combination of a force gauge and an accelerometer in one housing. 
This precision device allowed simultaneous measurement of 
acceleration and force very close to the driving point of the assembly. 
The impedance head also secured the threaded end of the 
compression bolt, providing a static pre-stress on the projector stack. 
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The opposite end of each stack was cemented to a solid brass slab five 
inches in diameter and one inch thick. 

2. Closed System Acceleration Calibration 

The following procedure was developed to match the 
displacements of the piezoelectric drivers, to ensure plug flow of the 
fluid in the chamber. The experimental set-up used is depicted in 
Fig. 8. An approximately constant voltage was simultaneously applied 
to each stack and the resulting accelerations were compared. This 
was done from 10 to 2000 Hz in ten-Hz frequency steps. At each 
frequency, the acceleration of one stack, as measured from the 
accelerometer in its impedance head, was adjusted to match the other 
stack's acceleration in magnitude (within +/- 2%) by manually varying 
the power amplifier. Phase adjustment was not required. A 
computer program (code listing appears in Appendix B) was written 
specifically to accomplish the tasks of: (1) measuring and comparing 
the acceleration magnitude of the two stacks at each frequency, (2) 
prompting the operator, when necessary, to adjust one stack to 
achieve a match within the desired tolerance, (3) recording the 
acceleration and force for each stack, and (4) stepping up the 
frequency to repeat the process. An HP 9000 Series 300 Computer 
commanded the HP 3421 A Data Aquisition/ Control Unit (DAC) and 
triggered the HP 4192A Low Frequency Impedance Analyzer. With a 
complete data base of force and acceleration magnitude and phase as a 
function of frequency, the displacement of the transducers was 
considered matched and known to within one percent. 
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Figure 8. Transducer Calibration Experiment Set-up. 



D. DIRECT METHOD EXPERIMENTS 
1. Introduction 

Several preliminary experiments were conducted using a fluid 
filled sample chamber with the objective to calibrate the direct 
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method apparatus. These experiments were performed with a single 
driver coupled to the test chamber in an "open system" configuration, 
as illustrated in Fig. 9. In Appendix A, it is shown that the radiation 
load presented to the open end is neglible in comparison to the 
moving mass load. Thus, no correction needs to be made to the 
theory in Chapter II for the open system. With the characteristics of 
the piezoelectric drivers established, a suitable means to measure the 
acoustic pressure resulting from the oscillating drive was developed. 
Finally, a set of calibration experiments was performed, the results of 
which were used to extract the moving mass of several "unknown" 
fluid loadings. 

2. Measurement of Acoustic Pressure 

Referring to Fig. 7, the original intention was to measure the 
acoustic force on the face of each driver using the impedance head. 
However, early experiments with various test masses showed that the 
moving mass could not be seen; i.e., the measured force was 
dominated by diaphragm stiffness, even though the diaphragm had 
been made as compliant as practical. To resolve this condition, a 
means to measure pressure in front of the piston, inside the sample 
cavity, had to be developed. Relocation of the pressure sensor 
permitted a stiffer diaphragm to be employed, and thereby reduce the 
effects of diaphragm flexure on the acoustic pressure in front of the 
piston. The initial diaphragms were made of 0.001 in. and 0.002 in. 
thick Mylar, and the final diaphragm was constructed from 0.003 in. 
thick beryllium copper. 
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Figure 9. Direct Method Equipment Set-up. 

Several alternative methods to measure the pressure inside the 
sample cavity were tested. Fig. 9 shows the equipment set-up. For 
these tests, the sample cavity was filled with various amounts of 
Fluorinert FC-72 non-polar electronic liquid (3M Corporation, St. Paul, 
Mn.). Fluorinert was used because it is non-conductive and it permits 
the selection of various viscosities that will be used in later (post 
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thesis) experimentation. First, a variety of Polyvinylidiene Fluoride 
(PVF 2 or PVDF) films, ranging from 0.001 in. to 0.01 in. thick 
(Pennwalt Corporation, Folsom, Ca.), were tried. Each sensor was 
bonded to the diaphragm. In every case, the resulting voltage 
response from the pressure sensor was contaminated by resonances of 
flexural waves in the diaphragm. Next, a similar series of trials was 
conducted using PVDF sensors suspended in the fluid column in the 
vicinity of the piston drive. Again, flexural waves dominated the 
response, in this case propagating in the sensors themselves. A 0.001 
in. thick electreted PTFE Teflon sensor was tested, with no 
improvement over the PVDF series. Finally, a small cylindrical tube 
piezoceramic hydrophone 6 mm in diameter was tested. The small 
hydrophone was selected because its lowest frequency resonance is 
well above the highest frequency of interest in this experiment. While 
the hydrophone lacked the sensitivity of the PVDF sensors, the 
measured pressure spectrum over the range 100 to 2000 Hz was 
substantially improved. 

Figure 10, raw data output from the HP 4194 for 80 ml added 
mass of Fluorinert, is an example of the voltage response of the 
accelerometer (lower trace) and the hydrophone (upper trace) as a 
function of frequency, plotted on logarithmic scales. The verticle 
scale for each trace spans five decades; the minimum and maximum 
values are listed. The verticle scale markings are for the 
accelerometer (A) trace. The frequency scale ranges from 30 to 3000 
Hz. The signal to noise ratio of the accelerometer output (lower trace) 
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Figure 10. Representative Response of Accelerometer ( "A: R" , lower trace ) 
and Hydrophone ( "B: T" , upper trace ) . 
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is very high, while the hydrophone output (upper trace) signal to noise 
ratio is poor at the lower frequencies. The accelerometer output is 
roughly proportional to f 2 for frequencies above 200 Hz, indicating a 
constant-displacement drive in this region. Below 200 Hz the 
acceleration is largely constant. For the most part, the hydrophone 
output signal is nearly proportional to the accelerometer output, the 
signature of a mass-like load presented to the driver. The peak in the 
hydrophone response at about 2000 Hz is a broad, contaminating 
resonance of unknown origin. As the use of a hydrophone is not 
compatible with the closed apparatus, it is regarded as an interim 
solution. 

Because of the low signal to noise ratio observed at low 
frequencies (10-100 Hz), considerable effort went into eliminating 
noise sources. Though a high driving voltage (^ 100 V) was desired, 
the available high voltage amplifier was found to introduce excessive 
noise into its output. A lower voltage output amplifier, Techron 7520 
Power Supply Amplifier (Crown International, Inc., Elkhart, In.), was 
selected for its low noise output. The Techron amplifier is limited to 
approximately 25 V rms output voltage. Pre-amplifiers were operated 
on battery power to avoid 60 Hz harmonic corruption in the pressure 
and acceleration signals. An exhaustive effort to eliminate ground 
loops among various equipment (in Fig. 9) proved to be most 
beneficial. 
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3. Drive Linearity 



Prior to collecting calibration data on the direct method 
apparatus, a characterization of the drive linearity was performed. 
Table 3.1 lists the accelerometer and force gauge outputs of the 
impedance head for applied voltages to the piezoceramic stack of 2V 
and 20V. The accelerometer and force gauge voltage responses were 
collected by the gain-phase analyzer. Data at each frequency were also 
evaluated for the presence of non-linearities using an HP 3580A 
Spectrum Analyzer. It was observed that the amplitudes of all 



TABLE 3.1. DRIVE LINEARITY 




Freq. ( Hz.) 


Acceleration ( pV) 




2 V 


20 V 


200 


16.5 


162 


498 


65.5 


651 


1005 


246 


2457 


2005 


1023 


10230 




Freq. ( Hz.) 


Force (mV) 




2 V 


20 V 


200 


5.94 


59.4 


498 


5.89 


58.9 


1005 


5.67 


56.7 


2005 


4.74 


47.4 



30 



harmonics were at least 30 dB below that of the fundamental driving 
frequency. As a result, it was concluded that the drive is linear. 

4. Calibration of the Direct Method Apparatus 

The objective of this final experiment was to obtain a set of 
calibration data for the apparatus which can be subsequently used to 
determine the unknown moving mass of a fluid contained in a rigid 
porous solid. An evaluation of the calibration data set was made 
possible by a comparison with results for known moving masses of 
fluid. For a nearly constant drive voltage of 25 V, calibration pressure 
and acceleration measurements were recorded as a function of 
frequency using the equipment in Fig. 9. Fluorine rt was added in 
increments of 10 ml , from 50 to 90 ml (the specific gravity of 
Fluorinert FC-72 is 1.68, the sound speed is 525 m/s). As with 
previous experiments using water, the viscous penetration depth of 
FC-72 is much smaller than the chamber sleeve (pore) radius; e.g., at 
100 Hz, S FC-72 = 100 M-. Plots were obtained for pressure and 
acceleration for the band 20 to 2000 Hz at each 10 ml increment for 
qualitative analysis, e.g., Fig. 10. Because a broad resonance around 
2000 Hz dominated the upper end of the frequency range, concurrent 
recordings of pressure and acceleration were taken at regular 
frequency intervals only over the range 100 to 1000 Hz. The data 
were collected using long time integration and 16 averages (features of 
the HP 4194), and the resulting calibration data set was stored in table 
format (a convenient function of the HP 4194). 
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Drive voltage uncertainty was less than 2% and added volume 
uncertainty was less than 4%. Most of the added mass error was 
attributed to evaporation of Fluorinert FC-72 at room temperature. 
Since extracted mass is a quotient of the two experimentally measured 
values, per Eq. (2.17), each known to better than 4%, overall 
experimental error is less than 6%. 

Using 60 and 80 ml data as the calibration control set, 
extraction of the moving mass for 50, 70, and 90 ml were 
accomplished at the selected frequencies from the following 
interpolation/extrapolation formula: 



Volume extracted ( co ) = 60 + (80 - 60) 



IE - JLE 






a'X 



'a '60 



— E .LIE 



L\V a /80 



'a'60J , 



(3.4) 



where, for a given mass loading and frequency: Vp is the output voltage 
of the pressure sensing hydrophone and V a is the output voltage of the 
accelerometer. The results of the extraction of the volume of moving 
mass, Fig. 11, are well within experimental error over most of the 
decade in frequency, 100 to 1000 Hz, and constitute a successful 
calibration of the direct method apparatus. Significantly, the ability to 
discern the volume of unknown added mass, or the resolution in ml, is 
5% or about 4 ml. 
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RESULTS : EXTRACTED MASS vs. FREQUENCY 




Figure 11. Calibration Results for the Direct Method. 
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IV. INTERPRETATION OF RESULTS 



A. COMPARISON OF TECHNIQUES 

It is important to compare the preliminaiy calibration results from 
the two methods and identify the more promising technique for future 
investigations with porous media of unknown complex mass. The 
impedance method is procedurally very easy to perform with the HP 
4194 Impedance/Gain-Phase Analyzer. However, because the 
extraction of moving mass is indirect, the impedance method is a 
technique which requires much finesse; a precise knowledge of the 
properties of the transducer is essential. What the direct method 
lacks in experimental ease is made up in its direct extraction of 
complex mass. 

Apparatus prototypes for both techniques were developed and 
employed successfully. The preliminary results of the impedance 
method, represented by Fig. 4, are comparable to those of the direct 
method, Fig. 11. Five per cent uncertainty in extracted mass over 
most of the frequency band 100 to 1000 Hz was established in both 
techniques. However, because the impedance method is plagued by a 
systematic deviation in extracted mass at higher frequencies, which 
increases in magnitude with increasing added mass, it is doubtful that 
moving masses greater than 50 g can be extracted using this 
technique with the present moving coil driver. The direct method 
exhibits no such error. Uncertainty in applying the direct method can 
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further be reduced by refinement of the design of the piston and 
surrounding diaphragm to reduce the effect of flexural waves. 



B. THE PREFERRED METHOD 

Based on the preliminary results of both experimented techniques 
to date, the use of the direct method for further investigation is 
recommended. While there is room for improvement in the apparatus 
of both techniques, the impedance method is limited by a driver that 
appears unsuitable for extraction of moving masses greater than 50 g. 
Conversely, it has been demonstrated that it is possible to extract a 
moving mass of approximately the value expected for fluid- filled 
porous solid samples (~ 100 g) with a resolution of several percent 
using the direct method. 
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V. CONCLUSIONS 



An experimental apparatus was developed by which the complex 
mass of a viscous fluid contained in a rigid porous solid can be directly 
extracted from acoustic pressure and acceleration measurements. 

Several preliminary experiments were conducted to verify the 
ability to extract the value of a known mass of fluid. The results for 
the extracted mass were accurate to within five per cent of known 
values over most of the frequency range 100 to 1000 Hz. 

A comparison of the preliminary results of the two techniques of 
obtaining complex mass demonstrated that the direct method was 
most suitable for use in further work investigating porous solid 
geometries. 
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APPENDIX A OPEN SYSTEM THEORY 



Calibration of both the direct and impedance methods was 
accomplished using the open system, single transducer/test chamber 
set-up of Fig. 9 (direct method) and Fig. 12 (impedance method). 
The effect of the radiation load presented to the open end on the load 
presented to the driver needs to be examined. Ignoring the radiation 
load for the moment, in the long wavelength limit, the impedance 
presented to the face of the driver by the liquid in the chamber alone 
is (Grant, 1988, pp. 8-10), 



where H is the mass of the liquid. Since the particle velocity is the 
same at the top and bottom of the liquid in the low frequency limit, 
the impedance seen by the driver is just the sum of the above, Eq. 
(A. 1), and that of the radiated sound wave in air: 



Substituting Eq. (A. 1) for Zraduo and the appropriate expression for 
Z,-,^ , assuming the worst case of ka = 1, into the right side of Eq. 
(A. 2) gives: 



^radpRivER — ^raduQ jtOM r 



(A. 1) 




(A. 2) 



^radoRrvER jwM r pcS ^ 



i DRIVER 



(A. 3) 



where S is the area of the sample at the fluid/air boundaiy. Since 



|pcS| 

IjcoMj 



« 1 



(A. 4) 



for the stated experimental conditions and dimensions, Eq. (A.3) may 
be approximated as: 



l^radi*rvEiJ — _ (A. 5) 

Thus, the theory in Chapter II does not need to be modified for the 
open system. 




Figure 12. Single Transducer /Test Chamber Open System. 
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APPENDIX B TRANDUCER CALIBRATION PROGRAM 
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DIRECT MEASUREMENT METHOD - COMPLEX MASS DENSITY 
CALIBRATION PROGRAM 












CONTROL 1 ,»Z;0 
PRINTER IS CRT 
PRINT USING " 
PRINT USING " //" 



i 



DorU'eys uiv 



written by R , A . MIRIDK , lest update Z 0 NO v 85 



01 Sr "Ensure Power ON to Analyzer, DAC , and all ancillary equipment 
WAIT E. 

PRINT " H 

01 Sr "Press CONTINUE when ready , , , " 

PAUSE 
DI5P " " 

• INITIAL CONDITIONS 

l*f»ir'»*»irr*t")Mr*r*f'frttif*r*»rT*'»»iMMr»tr'»*r*TrT’rinMt'riMrTifT»TrifT»TiMCTr 

OPTION BASE i 
INTEGER I , N p t s 

DIM Freq_nz( ZO0 ) 7 hccei7mf s ( ZOO ) , Acce 17aeg< ZOO ) ,ncce iomks ZOO % ' 

DIM ncce ioriegv ZOO ) ,ForceSmk s ( ZOO > 7 Force3deg( ZOO / 7 Porce 1 Gni . 5 ( 200 ‘ 

DIM Force 1 Odeg ' ZOO ) ,Ref_indbv< 200 ) ,Bm i nus_a7magdb ( 200 / 

DIM Bni nus_.a6me.gdb ( ZOO > , Acce i_7magv( ZOO ) ,Acce i_8megv' 200 > 

DIM Bni nus_aSmagdb ( ZOO / , Bmi nus_a. \ Omgtib v ZOO ) .For ce_5cn.agvi ZOO ) 

DIM Force_ 1 Ongv( 200 ) 

MASS STORAGE IS • ,700 7 \ " 

DEG 

usc_ie^e 1= , 5 ! Volts rns 

ASSIGN ennslyzer TO 7 \ 7 } Analyzer 



ASSIGN CDac TO 70S 
CLEAR ^Analyzer 
CLEAR §Dac 
WA IT 5 , 

STATUS ^Analyzer ; Srq_5 t atby te 
IF Sr o_s t a t by t e=64 OR Srq_s tat byte >S4 THEN 
DI5F " INPUT ERROR..." 

6 0 T G (5(0 
END IF 



unu 



nr * 1 
H r 3 4 Z i A 



' error exists 



I TERMINATE 
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a t 



w 

550 

540 

550 

5S0 

57G 

580 

530 

r+ r% r> 

a uu 
8 1 0 
820 
830 
540 
550 
850 
570 
680 
5 30 
700 
7 t 0 
720 
730 
740 
750 
760 



CALIBRATION PROCEDURE 

irirT-riririri 






rir-irnririntiti 






OUTFUT ©Analyzer ; "GL" ;Gsc_level ; "EN* 
OUTPUT ©Analyzer ; ,, N2T2vl " { 

! 

OUTPUT ©Analyzer; "EX " 

WAIT 4. 

Np t 5*0 

FOR 1^(0 TO 2000 STEP 10 
Np t s =Np t s+ 1 

Fr eq_r hz c I { 1 000 . 

OUTPUT ©Analyzer 
I WAIT i . 



1 FR 11 ; Fr eq_kh: 



External trigger , measure in dBU , 

Average on 



rrequencv in nz 



EN" ! Spot frequency set to I H 
! Allow settling time for impedance follower 
( operational amplifier 

BEEP 301 t . , » 5 

0 ISP "FREQUENCY IS ,1." Hz" 

WA IT 2 . 

OUTPUT ©Analyzer: "A6N2U1 " 

OUTPUT ©Anal v zer ; " EX " 

WAIT I . 

ENTER ©Ana Ivzer ; Re f_indb v( Np t s > 

OISF "Reference input is ,Ref_indb v( Np t s / , " dBU" 

WAIT 4 . 



Function. 1 measure REF input 

Trigger command 

Tate a ova ni age of AvERAGE- ing 



(01 Sr "Observe input on Oscop<( 

1 A5N2v 1 



770 


OUTPUT 


©Ana 1 y zer ; " i 


780 


OUTPUT 


©Oac ; " OrN " 


730 


OUTPUT 


©Dec : " CLS7 " 


800 


WAIT , 


0 1 


8 1 0 


OUTPUT 


©Ana ly:er; " i 


820 


WAIT 1 


* 


630 


ENTER 


©Ana 1 y zer : 5m 


840 


OUTPUT 


©Oac ; " 0rN7 " 


850 


OUTPUT 


©uac : " CLSS " 


680 


WAIT . 


0 1 


870 


OUTPUT 


©Ana 1 y zer r " ! 


88G 


WAIT 1 




8 30 


ENTER 


©Ana I y zer ; Bm 


300 


OUTPUT 


©uac ; " GrN6 " 


3 1 0 


i - - - 


dev to vo 


320 


Acce i_ 


7magv( Np t s )- 


330 


Acce 1_ 


6magv( Np t s /= 


340 


! - - - 


Volts to ( 


350 


Acce 17 


mf 5 ( Np t s )=Aci 


350 


Acceiomks ( Np t s >=Ac< 


370 


• - - - NO 


TE * THIS Ci 


380 


Ca l_d i 


f f -Acce 17mk s 


3S0 


OISF " 


Ac c < 



Ensure signal is not clipped 
Function* measure TEST input 
I Open all relays 

» Ensure relav closure 






in 4 32S7 S 



,ua i_o ir r 



( uuu 



40 



ei G) 



C3 CD (5 CD CD CD CD CD CD CD CD CD CD CD CD (D CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD C5 

in co co cn cd c'j fn in to n- co cn cd rJ 10 in to co rn cd cj rn nr in 10 t'- ct 

cd cd cd CD cd cd — — — — — — — — t'J rJ r J rJ rJ r j r j c j rj r j tn rn to \n tn ro rn rn K 



i <a t <a 
020 
030 



I 

WAIT 2 . 

if dirrvi . ( un oa I__d i f f <. . 3 THEN [ nccur-acy of 

DISF " WARNING.,. Acceleration ratio ls " ,Cai_aLff 
WAIT 4. 

INPUT ”uo you wish to try again? Y/N ” ,Answ3 
IF An5w3 s3,, Y" OR Answ3="y" THEN GOTO 730 
END IF 

OUTPUT ©Dac : " 0FN8 " 

OUTPUT ©Gac;"CLS3" 

WAIT .01 

OUTPUT ©Ana I y zer ; ” EX *' 

WAIT \ . 

ENTER sAna Lyzer ; an inus_a3magdb ( Np t s ) ,Force3deg( Np t s ) 
OUTPUT ©Qac;"0rN3" 

OUTPUT ©Dac;"CL3i0 u 
WAIT .01 

OUTPUT ©Analyzer EX " 

WAIT i . 

ENTER ©Analyzer t Sn mu5j i 0ngdb ( Np t s ) , Force i 0deg( Np t s ) 
OUTPUT ©Gac ; " OrN I 0 " 

[ - - - dSU to Uolts - - - I 

Force_3magv( Np t s )® I 0 . A ( ( Sninus_a3nagdb ( Np t s )+Re f_indb v( Np 
Force^i 0ngvt Np t 5 /* 1 0 . A ( ( Bn l nus_a \ 0ngdb ( Np t s ) + Re f_mdb v( Np 
5 ~ ~ volts to Newtons - - I 

Fore sank 5 ( Np t s )*Force_5nagv( Np t s >-*2,525 
Fore e t 0nks ( Np t s ) =For c e_ 1 0ngv( Np t s ) ■* 2 . 533 
?-***- NOTE J THIS CONVERSION APPLIES ONLY TO IH4 9237 & 33 
Fr e q_h z ( Np t s ) = I 
N E X T I 

STATUS ©Analyzer ; Srq_s tatbvte 

IF Srq_5 t atbyt e*S4 OR Srq_s tatav t e>S4 THEN * Error exists 

u I Sr "INPUT ERROR..." 

GOTO iSI® 1 TERMINATE 

END IF 



proceaur e 



t 5 / / ( ( 
t S ) ) / • 



4 1 



6) O 



I 400 
(4(0 
1 420 
1 430 
1 440 
(450 
\ 480 
(470 
1 480 
t 430 
1 500 
(510 
i 520 
1530 
( 540 
i 550 
1 55 0 
(570 
1 580 
1 530 
( 500 
1510 
1520 
(530 
(840 
I 550 
( 550 
( 570 
1580 
1 5 30 



I OUTFUT 



INPUT "Enter the OUTPUT file name. . . " ,Name$ 
CREATE BOAT Name* ,1300,8 
ASSIGN yFiie_name TO Nane$ 

FOR I = i TO Npts 

OUTPUT OF i le_name ; Freq_hz( I ) 

OUTPUT OF l ie^name ? Acce I7mks ( I / 

OUTPUT OF i ie_name ; Accel7deg( I ) 

OUTPUT OF L ie_name ; Acce 13mks ( I ) 

OUTPUT OF i Ie_name ; Acce I5deg( I ) 

OUTPUT OF l ie_name ; ForceSmks ( I > 

OUTPUT OFl I e_name ; ForceSdegv I ) 

OUTPUT OF l ie_name ; Force t Omke ( I ) 

OUTPUT OF l I e_name ; Force 1 0deg( I ) 

NEXT I 

ASSIGN OF 1 1 e_name TO * 






TERMINATE 






r -r -w i 






p-jf-r-* ibit-r-ie-r-f-ir-X’ritiff-xir + it-lt 






PRINT USING "0" 

PRINT US I NG " r / / " 

PRINT "Program conclusion. , . 
PRINT 

BEEP 405 , . 5 
BEEP 551 , . 5 
SEEP 835 5 

BEEP 132,1 * 0 
END 
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